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Abstract 

We evaluate the cohesive energies E& of four systems in which par- 
ticles move on a cylindrical surface, at fixed distance R from the axis. 
We find quite nonuniversal dependences of on R. For the Coulomb 
binding problem, E& is a monotonically decreasing function of R. For 
three problems involving Lennard-Jones interactions, the behavior is 
nonmonotonic; E& is larger at R = oo than at R=0; the maximum 
binding corresponds to R ~ 0.7a (the hard core parameter). Conse- 
quences of the enhanced binding are discussed. 

The discovery of carbon nanotubes has stimulated a rapid evolution of 
ideas, experiments and understanding concerning states of matter confined 
to the proximity of a cylindrical surface [jl], 0, [|. Examples of such systems 
include electrons present within a nanotube and atoms or molecules moving 
just outside or within such tubes. This paper reports unexpected behavior we 
have found in studies of four such systems: a +/— pair of charges bound by 
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the Coulomb interaction, a pair of atoms interacting with a Lennard- Jones 
(LJ) interaction, an ensemble of 4 He atoms, which condenses, and a low 
density fluid consisting of classical atoms. We assume that all particles move 
on a cylindrical surface, of radius R and infinite length. The assumption of 
surface confinement simplifies the calculations without sacrificing the basic 
physics. 

For each of these four systems, considerable attention has been directed 
previously to the investigation of two extreme limits of the present problem. 
The limit R = oo, here called "flatland", is that of particles moving on 
a plane, i.e., a two-dimensional (2d) problem. This has been extensively 
pursued in connection with both the 2d electron gas and monolayer films 
|fj, [5|, The opposite limit, R approaching zero, is here called "lineland", a 
Id limit. Matter in lineland has been explored for many years as an abstract 
problem |7J and has recently received particular attention in connection with 
the possible realization of Id phases within interstitial channels or external 
surface grooves on nanotube bundles 0. ||. |9], [K| 11], [12], |13|, [L4j [15|, 0. A 
logical question addressed in this paper is whether the properties of matter 
in " cylinderland" evolves smoothly (or even monotonically) between these 
limits as the value of R is varied. We find that the answer is "yes" in just 
one of the four cases and that intriguing behavior arises in all four cases. 

The first problem we address is the ground state binding energy E;,(R) of 
the Coulomb interaction problem. Consider a charge + e fixed in position 
(at z=0=O on the cylindrical surface) and a negative charge — e, mass m, 
which is free to move on the cylinder in the presence of an interaction — e 2 /r, 
where r is the charges' separation; the electron's position is at cylindrical co- 



ordinates (z, 0) and r = J z 2 + [2i?sin(</>/2)] 2 . The Schrodinger equation for 
the hydrogen atom is exactly solvable in 2d, resulting in a ground state wave 
function exp(— 2r/ao), with binding energy Eb(oo) = 4 Rydbergs and mean 
separation a /2 between the interacting charges. Here, ao is the Bohr radius. 
In the opposite limit of lineland, we encounter a well-known logarithmic di- 
vergence of Efe as R approaches zero; the particle "falls to the center of the 
attractive force" in Id. We expected the behavior of the cylindrical problem 
to shift between these limits when R becomes of order a . Figure 1 depicts 
the result of a three-parameter variational solution of the problem; details of 
the trial wave function will be provided in a future publication (as will those 
of the other calculations described here). One observes in Fig. 1 that E& 
increases monotonically as the reduced curvature C = a$/R is increased; the 
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known Id and 2d limits are accurately reproduced by this calculation. E& is 
not significantly affected by the nonzero curvature for C < f ; its value is only 
~ 1% higher at C — 1 than at C — 0. This weak dependence becomes com- 
prehensible when it is compared with the result from a perturbation theory, 
described elsewhere, where the curvature C is the perturbation. The resulting 
(lowest order) dependence on C is of the form E b (C) — E b (0) « (1/128)C 2 . 
As seen in the inset to Fig. 1, this function is consistent with the variational 
result, up to C = 1, explaining the initially weak dependence on C . 

The second problem we consider is the ground state cohesive energy (per 
atom) of a fluid consisting of 4 He atoms whose nuclei are confined to the 
cylindrical surface. The analogous problem has been extensively investigated 
in 2d |17| , for which it is relevant to superfluid films, and Id, where it has 
been studied for potential application to an interacting interstitial fluid |8], 
|9], [TI], O]. The cylinderland problem of 4 He has received some attention 



in connection with endohedral adsorption within nanotubes[lS]. In 3d, the 



cohesive energy of 4 He is E b = 7.2 K. In 2d, it is 0.85 K and in Id it is ~ 3 
mK. What R dependence is expected for E& on a cylinder? 

We have studied this problem variationally, using a Jastrow trial wave 
function, i.e. a symmetrized product of two-particle functions that prevent 
hard-core overlap of the atoms. The interatomic potential assumed in the 
calculation is the modern "Aziz" potential flT9| |. The results of this liquid 
state calculation appear in Figure 2. The dependence of E b on curvature is 
not monotonic. Indeed, the binding energy is a factor 3.7 higher near R=1.8 
A than in flatland. The nonmonotonic dependence on C and lowest value 
at C = oo stand in stark contrast to the monotonically increasing variation 
with C found for the Coulomb problem. 

The origin of the maximum binding energy near R=2 A is the minimum 
in the He-He interaction near r m ,„ = 2.9 A . The He fluid has a strong; 
binding if the geometry encourages the particles to have such a spacing. 
This is the case for the cylindrical geometry, as indicated by the following 
argument. We introduce a function called the specific area function a(r), 
defined as the area on the cylinder's surface at distance r from a specified 
point on the surface, per unit distance from this point. Letting this point be 
the origin, we have a(r) = J d 2 r' 5(\r' — r\). In the flatland limit, a(r) = 2nr, 
while in the lineland limit (R— > 0), a = 27rR. In the cylindrical case, a(r) is 
proportional to an elliptic integral, which exhibits a logarithmic divergence 
at r = 2R. The origin of this divergence is that a particle on one side of the 
cylinder has a divergent specific area at distance 2R; the sphere of this radius 
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is tangent to the cylinder. As a consequence, the system's energy is lowered 
when this distance is such that the interaction is strongly attractive. The 
optimal binding does not occur precisely when the minimum in the He-He 
potential coincides with the diameter of the cylinder, but instead at ~ 20% 
higher value of R. The difference arises from the zero-point energy of the 
system, which expands the nearest neighbor distance beyond r min . The same 
behavior occurs in 2d and 3d 4 He; the "nearest-neighbor" peak in the radial 



distribution function occurs at distance about 20% greater than r min pO . 

What consequences accompany the enhanced binding at this value of 
R? Typically (but not always) , strongly cohesive systems exhibit a relatively 
large speed of sound. Here, the sound propagation speed s is derived from the 
relation appropriate to longitudinal density fluctuations propagating parallel 



d 2 (E /N) 

to the axis of the cylinder: Ms 1 = p 2 — . Here, M is the atomic 

H dp 2 

mass and p is the Id density; the derivative is evaluated at the ground state 
density. According to our calculations, s = 230 m/s at the equilibrium 
density p = 0.3 A -1 for the optimally binding radius, R = 1.8 A . This value 
may be compared to the values s = 240,90 and 8.0(3.0) m/s in 3d, 2d and 
Id, respectively. We observe that the speed in cylinderland is significantly 
enhanced relative to both the Id and 2d values. 

Since the cylindrical fluid is a Id system, from the perspective of sta- 
tistical mechanics, it undergoes no phase transitions at finite T. There is, 
however, a T=0 transition as the system evolves from a liquid-vapor coex- 
isting ground state to a disordered fluid at nonzero T, with a singular heat 
capacity (proportional to 5(T)). The integrated specific heat is a monotonic 
function of the binding energy. This might be observable, due to inhomogene- 
ity in the system, as a smeared out maximum in the low T specific heat. The 
fluid's compressibility diverges as exp[£ , b /(/c B T)] at low T, which should be 
observable in an adsorption isotherm. However, we have thus far no definite 
calculations to compare with experiments. 

The third problem we address is the curvature-dependent binding energy 
Efc(R) of an atomic dimer on a cylinder of radius R. This is the two-body 
version of the liquid helium problem just discussed, except that here we treat 
the general case of atoms interacting with an arbitrary Lennard- Jones (LJ) 
interaction. By scaling distances relative to the hard-core diameter a and 

i i- ! „ . , /„ z ~ . r 

the energy relative to the well depth e [ z — — ; H — — ;r — — ; bb — — , 

V er a a e J 

the Schrodinger equation becomes: 
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Note that the ground state solution to this equation is determined by the 
boundary conditions and the value of the dimensionless de Boer quantum parameter 

h 2 

T] = — — - . A very large value of r\ implies a large zero-point energy and ab- 
sence of any bound state. Here, we focus on a specific question: what is 
the threshold value (rj = r] t ) separating those problems for which the dimer 
exists (77 < T) t ) from those for which it does not exist? The corresponding 
threshold value of rj t is known for the limiting cases of Id (j] t = 0.1788)0 
and 2d (r] t = 0.269) [2TJ. We have answered this question for cylinderland 
by computing the ground state energy variationally, identifying rj t from the 
point when the attractive potential energy becomes too weak for the dimer 
to be bound. Because the calculation is variational, the computed thresh- 
old rjt is a lower limit to the exact value. Our rjt results agree well with 
the known Id and 2d limits cited above. For the general cylindrical case, 
the threshold value is shown in Fig. 3. Qualitatively, it exhibits the same 
phenomenon as was seen in the 4 He liquid binding problem. That is, the 
binding is particularly large when the diameter of the cylinder is such that 
the interatomic interaction across the cylinder is strongly attractive. There 
is a small difference between the "optimal" radius, R op t, values for the two 
problems: R op4 /cx has the value 0.7 for the dimer problem and 0.77 for the 
4 He binding problem. 

Having in mind the interesting behavior of the threshold rj t in the cylin- 
derland, we next address the problem of the binding energy(E;,) of 4 He and 
3 He dimers on a cylinder. Employing the same variational approach as for rjt, 
we compute the ground state energy for these systems as a function of a/R 
( Fig. 4) . Not surprisingly, we find significantly enhanced binding energy for 
the 4 He and 3 He dimers at a particular range of values of R(~ 0.65 a). The 
enhanced binding in the cylindrical geometry is particularly dramatic for the 
3 He dimer, yielding an increase in E& of 7 orders in magnitude compared to 
the 2d limit for E&. 

Having established that the 3 He dimer is so strongly bound for radius R 
~ 0.65 a, we now discuss the ground state of the system of many 3 He atoms. 
No condensed liquid exists in either 2d or Id, but perhaps one exists in 
cylinderland. The existence of the stable dimer does not ensure the existence 
of a stable N-mer for any N > 2 (as was shown explicitly in the 2d case [p^l ) . 
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One scenario is that the ground state of the system is a gas of such dimers 
(analogous to H 2 at room temperature). A second possibility is that the 
dimers coalesce to form a liquid, analogous to liquid H2 between its triple 
point and its critical point. A third scenario is that the dimers dissolve into 
a many-body liquid ground state, analogous to liquid 3 He in 3d. None of the 
first two of these possibilities has been considered previously; the question of 
which phase, among the three candidates, is the actual ground state remains 
open. 

Finally, we address a fourth problem concerned with matter in cylinder- 
land — the second virial coefficient of a classical gas. By analogy with well 
known problems in other dimensions, we write a low density (high T) expan- 
sion of the Id pressure P = p 2 [ ^J- ) , where / is the free energy per particle, 

\dpj 

as P/(pk B T) « 1 + + ... , with 

B(T)= 1 -Jd 2 r[l-exp[-pV(r)}]. (2) 

Here, (5 = 1/(/cbT) and the integration is over the cylinder's surface. We 
assume the usual LJ form of interaction. The results(to be reported in detail 
elsewhere) exhibit similar qualitative behavior to that found for the preced- 
ing two problems, both of which involve an interaction with a well. As in 
the familiar 2d and 3d contexts, B is positive at high T due to the repulsive 
interactions and B < at low T, where attractive effects dominate. The 
Boyle temperature T^ is that for which B vanishes, meaning that (within 
this expansion) P for the interacting system is the same as that of a nonin- 
teracting gas at the same p and T. The results in Fig. 5 indicate that the 
Boyle temperature is the highest for a cylinder of radius R = 0.7 a. In the 
van der Waals theory of condensation, the critical temperature is 8T B /27. 
In that mean field theory, therefore, the critical temperature is highest for a 
cylinder of this "optimal" radius. While this transition does not occur in the 
exact theory, one expects the virial expansion to apply at low density and 
here the data of Fig. 5 imply an onset of the effects of attraction at higher 
T for the cylinder than for either Id or 2d limiting cases. Calculations using 
classical simulations for such gases should exhibit such an enhanced attrac- 
tion effect and extend the prediction to higher densities. We anticipate that 
crystallization of the classical gas will occur, with an oscillatory dependence 
of the energy on R (due to the size-dependent commensuration energy). Such 



6 



calculations are in progress in our group. 

In summary, we have explored both classical and quantum particles con- 
fined to a cylindrical surface. In the case of a Coulomb interaction, which 
lacks a preferred distance, the energetics of binding exhibits a monotonic 
trend as a function of curvature; there is an insensitivity of the cohesion to 
the value of the curvature, as long as it is small. In three problems involving 
interactions with a favored distance, quite distinct behavior was found. There 
is a particular range of values of R such that the binding is enhanced and 
this is explicable in terms of a phase space argument relevant to interactions 
which have well-defined potential energy minima. We believe that this dis- 
tinction is generic and should be applicable to other geometries. For example, 
matter confined to a spherical surface (small particles or pores) should have 
enhanced cohesion when the diameter is~ a. Other related problems merit 
investigation. One is the possibility of condensation of 3 He to a liquid. The 
factor of ~ 3 enhancement of the liquid binding energy shown for cylindrical 
4 He (relative to binding in flatland) suggests that there is a range of radius 
over which the lighter isotopic liquid should also bind. A suggestive argu- 
ment in support of that possibility is the fact that the 3 He dimer exists over 
an extended range of radii: _R>cr/2~1.3A , according to Fig. 3 (rj = 0.24 
for 3 He). For bose systems in Id and 2d, the dimer threshold coincides with 
the binding threshold for the many-body bound state liquid. While we do 
not know the corresponding criterion for self-binding in cylinderland; it is 
possible that the cylindrical liquid 3 He also exists. [] This would be a re- 
markable system to investigate since it exemplifies a novel Luttinger liquid. 
A primary application of these results is to carbon nanotubes, whose radius 
can be as small as 3 A . Indeed, the so-called cylindrical shell phase of He 
and H 2 corresponds to adsorption at radial distance 5r ~ 3 A inward from 
the carbon atoms [I8|. Thus, the pronounced effects found here for these 



gases at R ~ 2 A ( Figs. 2, 4 ) correspond to cylindrical phases within 
nanotubes of radius ~ 5 A . Experimental and other theoretical study of 
this size tube is worth pursuing. 
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1 Note that out-of-plane motion will enhance the likelihood of self-binding, as was found 
by Whitlock and Gordillo etaL@ for 4 He in 2d and Id and by Brami efaZ.Q 

for 3 He in 2d. The rms displacements are 0.25 to 0.4 A in these cases. 
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Figure captions 

Fig. 1 : Ground state energy of a ( H — ) Coulomb pair of charges, 
constrained on a cylinder as a function of a curvature parameter ao/R. The 
solid line is the variational result. The inset compares the variational and 
perturbation theory(PT) results. The energy is in Hartree units. 

Fig. 2 : Ground state cohesive energy (per atom) of a 4 He fluid on a 
cylinder. 

Fig. 3 : Threshold value of the de Boer parameter for existence of a 
dimer. 

Fig. 4 : Ground state energy of 4 He(circles) and 3 He(stars) dimers on a 
cylinder. 

Fig. 5 : Reduced Boyle temperature T* = /cgT/e . 
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